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ABSTRACT 



^ ", It is generally argued that present cosmological observations support the accelerating 

models of the universe, as driven by a cosmological constant or 'dark energy'. We 
| argue here that an alternative model of the universe is possible, which explains the 

current observations of the universe. We demonstrate this with a reinterpretation of 
the magnitude-redshift relation for Type la supernovae, since this was the test that gave 
a spurt to the current trend in favour of the cosmological constant. 



Introduction 



> 

There are only two instances in modern cosmology in which a genuine theoretical prediction has 
been directly confirmed by observation. The first of these was the theoretical prediction obtained 
directly from the solution of Einstein's equations by Friedmann (1922) and later Lemaitre (1927) 
that the universe is expanding (Hubble 1929). The second was the theoretical prediction made in 
the framework of the original steady-state cosmology (Bondi and Gold 1948; Hoyle 1948) that the 
universe should be accelerating. This has been demonstrated in modern times by Perlmutter et al 
; (1999) and by Riess et al (1998, 2000). 

Unfortunately while the latter was a clear cut theoretical prediction, the observers almost 
unanimously have chosen to interpret their detection of acceleration by maintaining the (standard) 
^ . and highly popular big bang model (SBBC) explaining the acceleration in terms of a non-zero 

cosmological constant. Because of this lack of balance in the discussion and the interpretation of 
the acceleration, we finally published a paper (Banerjee et al 2000) showing that the observational 
data available up to that time can very reasonably be interpreted in terms of the quasi-steady state 
cosmological (QSSC) model. In this model we had shown that the QSSC model provides a good fit 
to the above data provided one takes account of the absorption of light by intergalactic dust in the 
form of metallic whiskers. Dust of this kind is required by the QSSC model to explain the cosmic 
microwave background (Hoyle, Burbidge and Narlikar 1994a, 2000). 

Since that paper was published, new data have been obtained (Riess et al 2001) involving a 
type la supernova at a much higher redshift. 
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In order to give some level of balance to the claim that is now being associated with these 
discoveries (cf. Flippenko 2001), in this paper we extend our earlier discussion to see how well (or 
how badly) the quasi-steady state model fits the new data. However, we first begin with a short 
discussion of why the cosmological constant currently enjoys so much popularity. 

2. The advantages of A 

The recent popularity of the cosmological constant with the proponents of standard cosmology 
rests on the following aspects : 

1. The constraints of 'age' are relaxed. By choosing a A suitably it is possible to increase the age 
of the universe. Thus the age of a standard flat A = model is ~ 6.5 x 10 9 /iq 1 yrs, for Hq = 
100 ho km s" 1 Mpc -1 , which is too low for accommodating globular cluster ages of ~ 11 — 13 x 10 9 
yrs, for h « 0.7. By way of comparison, the age of a flat model with A = \/2>Hq = 0.7 is 
~ 13.5 x 10 9 yrs, for ho = 0.7. At the same time it should be stressed that the choice ho = 0.7 is 
rather arbitrary. While there are those who strongly advocate the model with this value, there is 
compelling evidence from Saha et al. (2001) that ho < 0.6. 

2. The use of A allows the metric distance to any given redshift to be larger than the corresponding 
distance (for the same redshift) for the A = model. This was demonstrated very clearly by the 
magnitude-redshift relation for Type la supernovae. The excess faintness of supernovae (over the 
predictions of the A = Friedmann models) was interpreted as the proof of a non-zero, positive A. 

3. Indirect checks on A / are provided by the structure formation scenarios and the observed 
inhomogeneities of the cosmic microwave background. Here the choice of a non-zero A helps by 
expanding the parameter space and thereby relieving pressures on structure formation scenarios 
from the various observational constraints. For a review of these and other observational constraints 
see Bagla eta al (1996) and Narlikar and Padmanabhan (2001). 

4. If one insists on a flat model (— as is implied by most kinds of inflation—) , the observed baryonic 
density of matter is only a few percent of the total cosmological density. In that case the rest (> 
90%) of all matter has to be dark and probably non-baryonic. The requirement of locating and 
identifying so much of matter in this dark and esoteric form is not so severe if we assume that a 
substantial part of it is taken over by A and is called 'dark energy' or 'quintessence'. Although the 
problem of bridging the gap between the predicted and the observed still remains, the gap is now 
to be explained more easily by an 'intangible' dark energy than by a more tangible dark matter 
particle. 
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3. The disadvantages of A 

The gifts brought by the cosmological constant come accompanied by a few unwelcome effects 
too. We name a few below. 

1. With the popularity of the inflationary model, the cosmological constant found a new life with 
its claimed origin in the force generated by transitions of the vacuum. It is this force that is 
supposed to drive the universe in a transient de Sitter-like mode. It is tempting to suppose that 
the present A owes its origin to the inflationary phase. However, Weinberg (1989) has commented 
on the extremely low value of A as required by the present cosmological observations compared to 
the primordial scenarios that drove the universe into an inflationary mode. If the inflation took 
place at the GUT epoch, the present value of A is too low by a factor ~ 10~ 108 . If the inflation 
took place at the quantum gravity epoch, the above factor is lower still at ~ 10~ 120 . It is easier 
to suppose that after inflation was over, the cosmological constant dropped to zero, rather than 
try to explain why we see its relic at such a low but highly fine-tuned value. One may recall that 
'fine-tuning' was a defect of the standard model, that the idea of inflation was invented to correct. 
To rationalize fine tuning, models with A oc H 2 , A oc S~ 2 or with A oc t~ 2 have been considered by 
Vishwakarma (2001) to allow for a steady drop of A from its high primordial value to the presently 
claimed low value. These models have a somewhat ad-hoc nature and certainly do not share the 
elegance of the overall structure of general relativity. 

2. Although it is tempting to place a larger and larger fraction of the total energy in the A - 
basket, beyond a certain level the exercise proves counter-productive. The stretching of space 
because of a large A makes gravitational lensing of large redshift objects, a much more frequent 
phenomenon than is observationally allowed. The presently favoured value of Ao = 0.7 is already 
bursting through this limit, and anything in excess may not be tolerated by the lensing data (Falco, 
Kochanek and Munoz, 1998) 

3. If one uses the 'best fit to observations' as a criterion for estimating A, then the angular size- 
redshift relation for ultra compact radio sources (Jackson and Dodgeson 1997) implies that A < 0. 
There is thus a contradiction between two cosmological tests of the universe going to large redshifts: 
the m-z relation gives A > whereas the 0-z relation requires A < for best-fit model. Given the 
uncertainties surrounding the cosmological 'standard candle' and 'standard rulers', we feel that the 
current enthusiasm for a positive cosmological constant may be premature. 

4. An alternative scenario : the QSSC 

If one is willing to move out of the framework of the standard cosmology, the quasi-steady state 
cosmology (QSSC) provides an already existing framework to understand the above cosmological 
issues. We briefly describe the salient features of this cosmology before coming to these specific 
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issues. 

The QSSC is a non-singular cosmology which arises from a Machian theory of gravity that 
gives field equations more general than general relativity (Hoyle, Burbidge and Narlikar 1995). The 
equations are: 

Rik - \g ik R + Xg ik = -^G[T^ + ^ reation ]- (i) 

Here we have taken the speed of light c = 1 and the first term on the right hand side is the usual 
energy momentum tensor of matter, which, in the case of "dust" (pressure p = 0) with density p 
and flow vector m, becomes 

rr tter = pu lUk , (2) 

while the second term on the right hand side denotes the contribution from a traceless scalar field 
C of negative energy and stresses with gradient C\ = dC / dx l : 

^creation = _ /[QC7fe _ ^QC 1 9ik }. (3) 

The coupling constant / is positive; but with its negative coupling to gravity or spacetime 
geometry, it generates repulsion and accelerated expansion of spacetime. We should point out here 
that the cosmological constant A in this theory arises from a Machian interaction and its sign is 
negative. 

For quantitative details of this model, see the paper cited above (Hoyle, Burbidge and Narlikar 
1995). General solutions of these equations in the homogeneous- isotropic case were obtained by 
Sachs, Narlikar and Hoyle (1996). The QSSC generally uses the simplest of these, which has a 
k = (flat) Robertson- Walker spacetime with a scale factor given by 

S(t) = exp(t/P) x [1 + n cos(27rr(t)/Q)]. (4) 

In this expression the function r(i) is very close to t, except near the maxima and minima of 
the scale factor. The parameters P and Q denote respectively the characteristic time scale of 
the exponential expansion and the time period of one cycle, while the parameter r\ satisfies the 
inequality 

< rj < 1. (5) 

Evidently for P » Q, the scale factor describes a universe with a long term exponential 
expansion together with short term oscillations in which S(t) never reaches the value zero. 
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The details of this cosmology are found in various papers by Hoyle, Burbidge and Narlikar 
(1993, 1994 a,b, 1995) as well as in their book {op.cit., 2000). The long term expansion in this 
model is accelerated, as in the original steady state theory. The repulsive force generated by the 
C— field is responsible for this feature. However, the C— field plays the crucial role of creating new 
matter without violating any conservation law. Although the scale factor describes a smoothed 
out universe, in actuality the creation is in the form of mini-explosions triggered by the strong 
gravitational fields of highly collapsed (black- hole like) objects. It is claimed that the creation 
phenomenon occurs at various mass scales, ranging from the stellar to the supercluster scales. 

Within the framework of the QSSC, the cosmic microwave background is explained as the 
thermalized relic starlight from previous cycles. The thermal temperature can indeed be worked 
out in terms of the known stellar activity and its present value can be derived as ~ 2.7K. How is 
thermalization achieved? The proposed mechanism is via metallic whiskers (Narlikar, et al 1997) 
which form a component of intergalactic dust. The formation and extinction properties of these 
whiskers have been discussed in the above paper. Typically the whiskers form from the metallic 
vapours ejected from supernovae where the metals are synthesized. Experiments have shown that 
metallic vapours on cooling, condense into elongated whiskers of ~ 0.5 — 1 mm length and ~ 10 _6 cm 
cross-sectional radius (Donn and Sears, 1963; Nabarro and Jackson, 1958). It is estimated that 
a dust density of ~ 10~ 34 g cm~ 3 is adequate to achieve the above effect. This model can also 
explain the observed anisotropy of the cosmic microwave background, including the peak at I ~ 200 
normally referred to as the Doppler peak in standard cosmology (Narlikar, et al 2002). Indeed it 
is this type of dust that helps explain the magnitude-redshift relation for extragalactic Type la 
supernovae. 



5. The magnitude-redshift relation in the QSSC 

From the above discussion, it is evident that three new aspects are introduced by the QSSC 
vis-a-vis standard cosmology when interpreting the m-z relation. There should be : 

(i) A negative cosmological constant, 

(ii) A repulsive force arising from a negative energy C-field, 

(iii) Cosmic dust which extinguishes radiation travelling over long distances. 

Of these (i) tends to decelerate the universe, while (ii) tends to accelerate it. The effect of (iii), so 
far as the Type la supernova m-z test is concerned, is to make the distant supernovae progressively 
dimmer. The effective magnitude-redshift relation is given by 



m(z) = M + 5 log(H d L ) + Am(z), (6) 
where M = M — 51ogi?o + 25, M is the absolute luminosity of the source and the luminosity 
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distance d^ is measured in Mpc. Here we have used the notation and definitions of Perlmutter, et 
al (1999). The quantity M has been referred to by these authors as 'Hubble constant free' absolute 
magnitude. The second term on the right hand side of this equation is the standard luminosity 
distance {di) term which can be computed for the QSSC. The product H^di is independent of 
Hubble's constant. The last term denotes the dimming produced by the intergalactic dust and may 
be explicitly written as: 



where, Ao = \/3Hq, fio = 8itGpo/3Hq and Q c fl = 8ttGp C: o/3Hq are the three parameters of 
this cosmology corresponding to the cosmological constant, matter density and the C— field energy 
density. Hq = 100 ho km s _1 Mpc -1 is the present value of Hubble's constant. The density p g o 
is the present intergalactic dust grain density in the form of metallic whiskers and k is the mass 
absorption coefficient. 

Chitre and Narlikar (1976) were the first to discuss the role of intergalactic dust in the m-z 
relation. At the time the general belief about the intergalactic medium did not include dust as 
a significant component and so the suggestion in the above paper that presence of even a small 
dust-component would enhance the apparent magnitudes was largely ignored. Recently Aquirre 
(1999) considered the role of dust in the m-z relation, followed by Banerjee, et al (2000). 

A comment is necessary on Eq (7). In the paper by Banerjee, et al (2000) this formula had 
appeared with a typographical error, with the factor y 2 missing in the numerator. The factor was 
nevertheless included in the actual calculation. The calculation redone by Vishwakarma (2002) 
has slightly lowered the value of x 2 P er degree of freedom from ~1.06 to ~1.00, thus marginally 
improving upon the goodness of fit claimed originally by Banerjee, et al (2000). 

Here we use this formula to compute the m-z curve as predicted by the QSSC. As reported 
by Vishwakarma (2002), the calculation is done in two different ways. The first is to use the low- 
redshift measurements to determine the two free parameters in equation (6), viz. M. and npgoh^ 1 , 
and then use these values to calculate \ 2 from the high-redshift measurements including SN 1997ff, 
the highest redshift supernova observed so far. Note that the model considered here is completely 
determined by Ao = —0.358 and rj = 0.811 as considered by Banerjee, et al (2000). The second 
method is to use the supernova 1997ff along with the earlier 60 supernovae, taking the high- and 
low-redshift points together, to get a best fit curve and then judge the goodness of fit. If the fit is 
relatively good, one can argue for the robustness of the theoretical relation. 

a) Method 1: Using the earlier work of Banerjee, et al (2000) we carry out the calculation of 
magnitudes at specified redshifts, using equation (6). Six of the sixty supernovae were dropped 
from the calculation as the deviations of the observed magnitudes from the expected ones are 
somewhat excessive. We will comment on this deviation within the framework of the QSSC later. 
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For the present we may mention that data on similar six supernovae had been omitted in obtaining 
the best-fit models of standard cosmology by Perlmutter, et al (1999) who had also considered 54 
supernovae in their primary fit C. 

With the constraint p g o > 0, the sample of 16 low redshift supernovae gives Ad ~ 24 with a 
wide range < Kpgoh^ 1 < 10 as the best-fitting solution, where x 2 does not change significantly. 
For example, A4 = 24.03 with Kpgoh^ 1 = 6.5 gives a value x 2 P er degree of freedom of 0.8 which 
represents an excellent fit. To fix ideas, A4 = 24 gives the corresponding absolute magnitude at 
peak luminosity of the supernova as M = —19.49 for ho = 0.6 and M = —19.16 for ho = 0.7. These 
values are in the right region for SN peak magnitudes (Saha, et al 2001). When we add the new 
point SN 1997ff (z = 1.755 ±0.05, m = 25.68 ±0.34) to the older list of 38 high redshift supernovae 
with z < 0.8, the new sample gives a x 2 P er degree of freedom of 1.49 which, though larger than 
before, is still acceptable as a reasonable fit. It would be worthwhile to mention that the supernova 
1997ff is claimed to be on the brighter side due to gravitational lensing effect (Moertell et al, 2000). 
If this is the case then an increase in the measured apparent magnitude of this supernova by only 
as much as 0.5 improves the fit considerably by giving x 2 per degree of freedom of 1.22. Likewise, 
a change in the redshift of this farthest supernova from 1.755 to 1.65 improves the x 2 P er degree 
of freedom to 1.33. We demonstrate this to underscore the need for caution in interpreting data 
in terms of models till (a) we have data on more supernovae available and (b) we are more certain 
about their magnitudes and redshifts. 

b) Method 2 : Here we include the supernovae 1997ff along with the 54 supernovae from earlier 
samples in obtaining the best-fit curve. Naturally the fit of such a curve to the data will be better 
than that in the first case. Here we report an exercise carried out by one of us (Vishwakarma, 2002) 
in which the models are fitted in the QSSC parameter space of > Ao > —0.32 and p g o in the 
range 3 — 5 x 10 _34 g cm" 3 . The region enclosed by the tilted U-shaped curve in Figure 1, has the 
QSSC models with parameters that yield values of x 2 P er degree of freedom in the range 1.16 to 
1.26, that is, values indicating a satisfactory fit. Even the model chosen by Banerjee et al (2000), 
viz. Ao = —0.358 yields a reasonable fit with the x 2 value per degree of freedom of 1.27. Moreover, 
if we assume that the supernova 1997ff is observed on the brighter side, as we have mentioned 
earlier, then an increase in its magnitude by 0.5 improves the x 2 P er degree of freedom to 1.12. 
Likewise, a change in its redshift from 1.755 to 1.65 reduces the x 2 P er degree of freedom to 1.19. 
In Figure 2, we have shown the fitting of this model (Ao = —0.358, rj = 0.811 with x 2 per degree 
of freedom = 1.27) to the observed data points and compared it with the favoured standard big 
bang model (f^totai = !)■ It should be noted that the Perlmutter et al' data taken together with SN 
1997ff, do not give the standard flat model as the best-fitting model: there the best-fitting model 
is obtained with Qo = 0.87, Ao = 1.51 (i.e., ^totai = 2.38) with a x 2 value per degree of freedom 
of 1.09. However, if we concentrate our attention on the flat model (f2 tota i = ^o + Ao = 1), which 
is currently favoured, the minimum x 2 is obtained for the model f^o = 0.34 with a x 2 value per 
degree of freedom of 1.17 (Vishwakarma, 2002). 
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Refering back to our opening remarks in this paper, it is interesting to compare the old steady 
state model with the classical Einstein deSitter model (0,q = 1, Ao = 0) for the data of Perlmutter, 
et al (1999). The respective value of x 2 per degree of freedom for the steady state model is 1.42, 
compared to the value 1.75 for the Einstein deSitter model. This relative performance may be 
contrasted with the situation prevailing three decades ago when the m-z data on galaxies were 
claimed to support a decelerating universe and were used to rule out the steady state model on the 
grounds that the model predicts an accelerating universe. 

6. Concluding remarks 

As was originally stated by Banerjee, et al (2000), the adjustable parameter in the above 
exercise is p g o, the density of intergalactic dust. By requiring the fit to be 'best' we determine 
the value of p g Q. Had the theory been on the wrong track, this value could have come out vastly 
different from the value needed to thermalize the microwave background. That it comes out in the 
correct range must be considered a point in favour of this approach. 

It will take us too far from the focus of the present paper to elaborate on the details of the 
microwave background in the QSSC. In a separate paper (Narlikar, et al 2002) devoted to the 
spectrum and fluctuations of the MBR in the QSSC we have shown that the MBR anisotropy 
from rich clusters can be revealed by a peak in the power spectrum around I = 200 while smaller 
peaks are revealed in the I ~ 600 — 900 range from groups of galaxies and perhaps weaker peaks at 
I > 3000 from individual galaxies or compact groups. 

In the QSSC, the whisker dust is produced and ejected from the vicinity of supernovae, as 
discussed by Narlikar et al (1997). Thus it is quite likely that dust near the parent supernova 
may show varying degree of concentration. This can lead to relatively large fluctuations in the 
supernova magnitudes, and may account for the large deviations for the six supernovae omitted 
from the study. 

In principle, despite such fluctuations, it is possible to use the m — z test to distinguish between 
the standard A-based approach and the above QSSC approach. As we go to higher redshifts, the 
effect of acceleration as implied by the cosmological constant will diminish, while that of the 
intergalactic dust will become greater. However, great care needs to be exercised before drawing 
any conclusions from such observations. As the example of the m — z test for galaxies showed in 
the 1970s, uncertain evolutionary effects tend to dominate at high redshifts. In addition, we have 
seen that even small changes in the measured values of m or z can greatly alter the conclusion. 
Finally, as we go to higher redshifts, gravitational lensing may selectively amplify the images and 
thus push the observed brightness to values higher than normal (Lewis and Ibata 2001). Hence it 
is prudent to await the accumulation of more type la supernovae in the high-z basket. 

Geoffrey Burbidge would like to thank one of his collaborators (JVN) for hospitality afforded 
him at IUCAA. 
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Fig. 1. — The tilted U-shaped contour shows a parameter-space of interest in the QSSC. The points 
of the contour have x 2 per degree of freedom in the range 1.16 — 1.26. The coupling constant k 
and the present density of the whisker grains p g o are measured, respectively, in units of 10 5 cm 2 
g _1 and 1CT 34 g cm~ 3 . [Adapted from Figure 6 of Vishwakarma (2002)] 
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Fig. 2. — Hubble diagram for 55 supernovae is shown. The theoretical curves represent the best- 
fitting flat QSSC (A = -0.358, rj = 0.811) and flat SBBC (fi = 0.34) models. The six omitted 
points are shown by open circles. 



